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One Sentence Summary: 

The dimensions and the age of Saturn’s main rings may be controlled by Yarkovsky forces, 
which also cause the particles in the C ring to fall on the planet, where the lu min ous fireballs 
create a Ring of Fire around Saturn’s equator. 
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Abstract 

The dimensions of Saturn’s A and B rings may be determined by the seasonal 
Yarkovsky effect and the Yarkovsky-Schach effect; the two effects confine the rings 
between -1.68 and -2.23 Saturn radii, in reasonable agreement with the observed values of 
1.525 and 2 . 261 . The C rin g may be sparsely populated because its particles are transients 
on their way to Saturn; the infall may create a luminous “Ring of Fire” around Saturn’s 
equator. The ring system may be young: in the past heat flow from Saturn’s interior much 
above its present value would not permit rings to exist. 
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What determines the dimensions of Saturn’s ring system? Are the rings young or 
old? Below it will be shown that the interplay of two Yarkovsky forces predict that the outer 
edge of the A ring and the inner edge of the C ring should be at certain distances from 
Saturn, and that these agree well with observation. The Yarkovsky forces also imply that the 
rings are relatively young. 

The interplay is between the seasonal Yarkovsky effect, which makes the orbits of 
ring particles decay, and the Yarkovsky-Schach effect which, for random spin axis 
orientations, always makes the orbits expand. The seasonal and diurnal Yarkovsky effects 
have been extensively treated elsewhere ( 1 - 4 ); a discussion of the Yarkovsky-Schach effect 
follows. 

The Yarkovsky-Schach effect is about how solar thermal forces on the ring particles 
behave when the particles enter and exit Saturn’s shadow (5,(5). The basic idea is 
qualitatively illustrated in Figure 1. The particle absorbs sunlight and then reemits it in the 
infrared. The infrared photons carry away more momentum from the hotter hemisphere 
than the colder one. The magnitude of the forces is indicated by the size of the thick arrows. 
When the particle dives into the shadow, the solar heating turns off, but the force takes some 
time to die away, due to the thermal inertia, and the particle continues to get a kick which 
increases the semimajor axis. As the particle emerges from the shadow, the force now acts 
against the motion so as to decrease the semimajor axis; but because the particle has cooled 
off, it takes some time for the particle to warm up and the force to turn back on. The net 
effect is a positive along-txack acceleration which increases the particle’s distance from 
Saturn. There are geometries for which the semimajor axis decreases; but for random spin 
axis orientations, the along-track forces do not average out to zero, resulting in a secular 
increase in the semimajor axis: the particle moves away from Saturn, increasing its angular 
momentum. A quantitative treatment of the Yarkovsky-Schach effect follows. 
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Figure 1 . Schematic of the Yarkovsky-Schach effect A ring particle orbits around Saturn; the view is 
looking up at Saturn’s south pole. The particle’s spin axis lies in the orbital plane and points to the Sun. 
The line across the particle is its equator. Sunlight comes in from the left and heats the sunward 
hemisphere, causing a thermal force (thick arrows). The size of the arrows indicates the magnitude of the 
force. The forces do not balance around the orbit because of the shadow and thermal inertia. The imbalance 
causes the orbit to expand secularly. There are also spin axis positions which cause the orbit to decay. 
However, averaging over all possible spin axis orientations causes the orbit to expand, as shown in the 
text 

The amount of sunlight falling on the particle will be given by b(<f>)F Sun , where F Sun = 
1371 W m' 2 x (1 AU/9.58 AU) 2 = 14.94 W in 2 is the solar insolation at Saturn’s distance 
(9.58 AU) from the Sun, and b(<p) is the “boxcar” function, which is equal to 1 when the 
particle is in full sunlight and 0 when it is in Saturn’s shadow. Also, <p ~ nt is the angle in 
the equatorial plane from the x-axis (see Figure 2), where t is time and n = {GMJa 3 ) 111 is the 
mean motion, with M s = 5.688 x 10 26 kg being the mass of Saturn and G = 6.67 x 10' 11 m 3 
kg" 1 s' 2 being the universal constant of gravitation. Simple Fourier analysis shows b(<p) to be 
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b($) = b 0 + 2 6,- cos U<P) 

7=1 

where £> 0 = 1 — A and bj = 2 sin[/(l — A)7r]/(/7r) for all other terms, with A being the fraction 
of the orbit which is shadowed. The parameter A ranges from 0 (outside the shadow) to 1/2 
(hugging the planet). 

Let x, y, and z be the unit vectors along the respective axes in the reference frame 
shown in Figure 2. The Sun is stationed in the x-z plane, so that the unit vector r s pointing 
from Saturn to the Sim is given by 

r s = sin0 s x + cos 6 S z (1) 

where 0 S is the Sun’s colatitude. Likewise, let s be the unit vector in the direction of the 
particle’s spin axis, so that 

s = (sint^ cos X p ) x + (sint^ sin X p ) y + (co s6 p ) z 

where 6 p is the axis’ colatitude and X p its east longitude. r s and s change slowly enough so 
that f s • s can be considered constant over one revolution of the particle about Saturn. 



X 


Figure 2. Coordinate system centered on Saturn. The z-axis points in the direction of the positive spin 
axis. The Sun lies in the x-z plane at a colatitude of 6 S . A ring particle P orbits in the equatorial plane at a 
distance a from Saturn’s center, and makes an angle <p=nt with the jc- axis, where n is the mean motion and 
t is time. 

The only term in b((f>) which will give a secular (i.e., long-term) evolution for the 
semimajor axis is the j = 1 term, because it depends on cos 0. The Y arkovsky-Schach 
acceleration f YS of a spherical particle will then be given approximately by 


f re = - (1 - A. v )B'b 1 F Sun cos (<p - 5)(r s • s) s 
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where B’ and 8 are an amplitude and a lag angle determined by the size and thermal 
properties of the particle and A v is its albedo in the visible. In effect the thermal inertia 
rotates the shadow by an angle 8. The along-track acceleration S YS is given by 

S ys = (f ra • t) , where t = - sin 0 x + cos 0 y is the unit vector tangent to the orbit. Working 
out all the geometrical factors yields and averaging 0 over one revolution about Saturn 
yields $YS being proportional to 

- 2 (r s • s) (s • t) cos (0 - 8) 

= + sin 6r sin 2 6„ cos 2 X„ sin <5 + cos#,, sin <9 cos6L cos X„ sin 8 

8 p p 8 p p p 

- sin 6 S sin 2 6 p sin X p cos X p cos 8 - cos 6 S sin 6 p cos 0 p sin X p cos <5 . 

For time scales longer than the orbital period, all but the first term on the right side 
of the above expression will drop out when averaging over longitude X p , because of 
precession. Note that the first term is greater than or equal to zero, so that S rs > 0; the 
Yarkovsky-Schach effect tends to increase the semimajor axis a of the particle. Thus the 
Yarkovsky-Schach effect pumps angular momentum into the rings. Averaging this term 
over all possible spin axis orientations gives the double average < S YS > = +(1-A V ) F Sm 
B'b^sm 6 S sin 813. Lastly, because Saturn’s equator is tilted by 26.73° with respect to its 
orbital plane, hjSin 6 S must be averaged over one Saturnian year to account for the changing 
shadow geometry. 

The diurnal Yarkovsky effect will average out to zero over all spin axis orientations. 
The averaged along-track acceleration due to the seasonal Yarkovsky effect is < S w > = — S 
( l-A m )B ’F Sun R 2 sin <5 /6a 2 for circular, equatorial orbits (3), where R s = 60,330 km is the 
radius of Saturn, a is the semimajor axis of the particle, A m is the albedo in the infrared, and 
the factor E = 1.79±0.09 is the amount that Saturn’s thermal radiation exceeds the solar 
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input (7). The subscript “TD” stands for “Thermal Drag,” the ori ginal name for the 
seasonal Yarkovsky effect (7). 

For spherical particles B’ is given by 


B' = 


9cpK 


cf + c\ 

c? 


in the theory which relates temperature linearly to flux, where c = 3 x 10 8 m s' 1 is the speed 
of light, o- 5.67 x 10' 8 W m' 2 K" 4 is the Stefan-Boltzmann constant, p is the particle’s 
density, K is its thermal conductivity, and T 0 is its average temperature, while £ m is its 
infrared emissivity. Also, tan 8= C/C p where C,, C 2 , and C 3 are functions which appear in 
the spherical Bessel function formulation of the temperature and depend in a complicated 
way on the particle’s radius R and thermal properties (5). 

The total along-track acceleration is < 5 > = < S YS > + < S ro >. Figure 3 plots < S 
> for a water ice “fluff ball” with the particular values R- 1 cm, T 0 = 90 K, p = 380 kg m 3 , 
£ m =1, A v - A m = 0, and K- 2.2 x 10 -4 W m _I K' 1 (8,9); but where the curve crosses the 
horizontal axis is independent of these parameters so long as A v = A^ because the 
Yarkovsky forces vary in the same way for both effects. The curve crosses the axis in two 
places, at distances a BC and a A0 . The region in between is where the net force on the particle 
is greater than zero and where rings would be expected. Beyond a A0 the seasonal 
Yarkovsky effect predominates and a particle would be forced back inward into the ring 
system. At distances less than a BC the seasonal effect again predominates, and particles 
would fall onto Saturn. 
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Figure 3. Averaged along-track acceleration < as a function of a. Semimajor axis a is measured in 
units of Saturn’s radius R s . The acceleration is the sum of the Yarkovsky-Schach effect and the seasonal 
Yarkovsky effect For Saturn’s current heat output < S> is positive between 1.68 R s and 2.23 R s . For 
values of a > 2.23 R s ring particles would be forced back into the rings, while those with a < 1.68 R s 
would fall on Saturn. 


For the above parameters, the Yarkovsky forces predict the outer boundary to be at 
a A0 = 2.23 R s \ the observed outer edge of the A ring is at 2.267 R s , where there is a 7:6 
resonance, and the two values agree to with 2%. The Yarkovsky inner boundary is at a BC = 
1 .68 R s ; the boundary of the B and C rings is at 1.525 R s (10,11), so that these two values 
agree to within 10%. Apart from the resonance, the differences might be attributable to 
unequal visible and infrared albedos or to E being less than 1.79; S = 1.70, for instance, 
yields a A0 - 2.34 R s and a BC = 1.54 R s . It is also of interest to note that the net acceleration 
reaches a maximum at 1.98 R s , which is close to the center of the Cassini gap at 1.987 R s . 

The graph given by Figure 3 indicates that the outer edge a A0 is the most stable 
position for the particles. So why haven’t all the particles migrated there? The answer may 
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be that the A and B rings axe “full”: collisions make the rings want to spread out into a 
near monolayer of particles which are as tightly packed as ring dy nami cs allows. The 
restoring forces at the outer edge of the A ring increase with distance until about 2.80 R s 
(where they start leveling off), so that spreading much beyond the outer edge of the A ring 
is difficult. On the other hand, at the inner boundary a BO once particles cross the inner edge 
of the B ring into the C ring, the seasonal .Yarkovsky effect drags them inward, and it is an 
easy trip to Saturn. Thus one can envision the following “reservoir” model for the A and 
B rings: they are filled with particles. Collisions cannibalize the larger particles and make 
smaller particles. The new particles elbow the older ones for room. The older ones find it 
easiest to spiral into the planet, rather than extend the A ring outward. 

This may also explain why the C ring, which is located where thermal drag 
predominates, is more sparsely populated than the A or B rings: the particles in the C ring 
are transients, on their way to Saturn. If this is so, then there might be a continuous rain of 
icy meteoroids all around the planet’s equator, whose l umin ous fireballs would give the 
impression of a “Ring of Fire.” The Ring of Fire might be visible to Cassini’s cameras 
when imaging the equator on the night side. 

Figure 4 graphs semimajor axis a vs. time t for particles with radii R- 0.1 cm, 1 cm, 
and 10 cm using the equation 


da _ 2(S) 
dt n 

(12). The particles start their journey at 1.70 R s . It is assumed that the only forces 
operating on the particles are Saturn’s gravity and the Yarkovsky forces. The orbit evolves 
the quickest for the 1 cm particle, which has nearly ideal properties for the operation of the 
Yarkovsky effects. The rings consist predominantly of particles in the size range 1-10 cm, 
close to the ideal. 
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TIME (10® YR) 

Figure 4. Semimajor axis a as a function of time t for spherical particles with radii 0.1 cm, 1 cm, and 


10 cm, assuming the particles are “fluff balls” and act under Yarkovsky forces and Saturn’s gravity only. 


No collisions or gravity from the moons are assumed. The particles are initially at 1.7 R s . 


Collisions will slow down the rate of Yarkovsky evolution by degrading the 
coherence between thermal forces and the orbit However, because both the Yarkovsky- 
Schach effect and the seasonal Yarkovsky effect have the same functional dependence on <p, 
they will decohere the orbit-spin coupling in lockstep, so that the positions of a BC and a A0 in 
Figure 3 will remain unchanged. 

Estimates of how often collisions take place range from 1 to 1000 times per orbit 
(13-15). Obviously bumping is not serious for sparse, and probably not top serious for 
close-packed particles. To see this for the close-packing, assume the particles collide 1000 
times per orbit. This is about once every 30 seconds. The particles most likely rotate with 
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angular speeds close to synchronous rotation (16-19). If each collision causes the spin axis 
to precess in a new direction at the same rate as rotation, the axis will rotate only through 

0. 4. before being struck again. With Brownian motion, the angular position will grow like 
N 1 ' 2 , where N is the number of collisions. So after one orbit and 1000 collisions, the axis 
will have moved through only 12°, which is not terribly significant for the direction of 
Yarkovsky forces. 

So probably somewhere in between sparse and close-packed is the regime where 
coherence degrades the most This leads to the idea of Yarkovsky-collision feedback: when 
the particles huddle together and collisions become frequent, Yarkovsky efficiency improves 
and spreads the particles out. Hence Yarkovsky may help explain the spacing between ring 
particles. Collisions, and the gravitational perturbations from Saturn’s many moons, await 
future research. 

Yarkovsky may constrain the age of Saturn’s rings. At present E = 1.79+0.09; but 
in the past E was larger, due to the higher heat flow inside Saturn (20), so that the seasonal 
Yarkovsky effect was greater and the zone where < S > is positive was smaller; a BC and a A0 
were closer together. At E = 1.90, a BC = a AO and the rings are unstable; the seasonal 
Yarkovsky effect predominates at all distances, and ring particles must fall on Saturn. Thus 
the Yarkovsky effects imply the rings are young, perhaps on the order of -10 s years old, in 
agreement with other mech anisms (15). Possibly they were created in the cometary shower 
which occurred 36 x 10 6 years ago (21,22). 
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